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A detailed analysis of the classical nonlinear dynamics of a single driven square potential barrier 
with harmonically oscillating position is performed. The system exhibits dynamical trapping which 
is associated with the existence of a stable island in phase space. Due to the unstable periodic orbits 
of the KAM-structure, the driven barrier is a chaotic scatterer and shows stickiness of scattering 
trajectories in the vicinity of the stable island. The transmission function of a suitably prepared 
ensemble yields results which are very similar to tunneling resonances in the quantum mechanical 
regime. However, the origin of these resonances is different in the classical regime. 



PACS numbers: 05.45.Ac, 05.45.Pq 



I. INTRODUCTION 

Harmonically driven potentials appear in many ar- 
eas of modern physics, particularly in mesoscopic elec- 
tronic semiconductor devices and other micro- and nano- 
structures driven by external voltages or applied laser 
fields. They also play a role for ultra-cold atomic wave 
packets exposed to optical barriers and photo-induced 
dynamics in strong laser fields or dissociation processes 
of molecules on solid surfaces. The strong external driv- 
ing of the system typically leads to nonlinear quantum 
effects and chaos in the corresponding classical systems. 
Two archetypical potentials have been investigated in de- 
tail in the literature, the driven potential well and the 
driven potential barrier. The periodic driving can be 
either a driving of the height or of the position of the 
potential. In the following, we will give a short overview 
of the known features of these systems. 

An early study of a vertically oscillating rectangular 
potential barrier, i.e. a potential with harmonically os- 
cillating height, in [lj aimed to derive an expression for 
the tunneling time through potential barriers. Particles 
interacting with a driven potential barrier can absorb or 
emit quanta offuv, where u> is the driving frequency. This 
leads to frequency-dependent resonances in the tunnel- 
ing probability through the vertically oscillating barrier, 
as shown in [2| for a rectangular barrier and a raised co- 
sine potential. These resonances can, according to Q , be 
interpreted as poles of the scattering amplitudes in the 
complex plane. If the potential barrier is delta-shaped, 
it possesses a set of leaky bound states, which have been 
detected in [ij by locating the complex energy poles of 
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the transmission amplitude. These semi-bound states 
lead to additional resonances in the transmission func- 
tion. Other works analyzed the tunneling through a ver- 
tically driven Gaussian-shaped barrier [5( and through 
the "Eckart" -barrier [6|. Both found an amplification of 
tunneling for intermediate frequencies, but did not de- 
tect a resonant behavior, which is due to the parameters 
chosen in these works, according to Q. A vertically os- 
cillating rectangular barrier, enclosed in a rigid poten- 
tial box, has been studied in the classical regime, see 
refs. HBO- The phase space of this system is mixed, 
with a purely chaotic layer at low energies, KAM-islands 
at intermediate energies and invariant spanning curves 
at high energies. The distribution of transversal times 
through the oscillating barrier is asymptotically algebraic 
with an exponent 7 = — 3. The Lyapunov exponent in 
the chaotic sea was also calculated for this system in 
as a function of the parameters. The Lyapunov exponent 
changes abruptly whenever an invariant spanning curve, 
separating different parts of the chaotic sea, is destroyed. 
If the driving is stochastic instead of harmonic, the sys- 
tem has no invariant spanning curve for high energies and 
the particles exhibit normal Fermi- acceleration ([10(). 

A laterally oscillating potential is often the result of a 
"Kramers - Henneberger" transformation of an ac-driven 
static potential. For high driving frequencies the tun- 
neling probability through an ac-driven rectangular bar- 
rier shows resonances at small energies for which the 
static barrier would be entirely intransparent, see (llj 
and [7[ . This can be described as resonant tunneling into 
quasi-stable bound states of the effective time-averaged 
potential, which has a double-barrier structure. For in- 
termediate frequencies, the scattering process is domi- 
nated by inelastic processes and strong sidebands in the 
energy spectrum Q- Such a resonant behavior was not 
found in the ac-driven Gaussian-shaped barrier Q and 
the "Eckart" -barrier [6]. Instead, these two systems ex- 
hibit phase-sensitive tunneling resonances for intermedi- 
ate frequencies, which can be explained by an increase 
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in the relative kinetic energy of incoming particles when 
the barrier approaches them. A moving potential barrier 
can also be used to tailor wave packets or to split an ini- 
tial pulse into several well-separated coherent pulses, see 

tutu. 

Potential wells with oscillating bottom have been stud- 
ied in 14] and In the classical regime, the os- 
cillating square well is pseudo-integrable and therefore 
not a chaotic scatterer (see also [la]), whereas the os- 
cillating smooth well has a stable KAM-island in phase 
space which leads to chaotic scattering. In the quan- 
tum mechanical regime both the smooth and the rectan- 
gular oscillating well accommodate quasi-bound Floquct 
states. Multi-photon processes couple incoming particles 
to these states, leading to resonances, which are visible 
as dips in the transmission function and as peaks in the 
dwclltimc. 

The tunneling through the laterally oscillating square 
well has resonances due to quasi-bound Floquet states as 
well ([3] and [IB])- However, due to the lateral driving 
these states belong to an effective double-well potential 
which leads to the formation of tunneling doublets. In 
the classical regime, the ac-driven square well has a sta- 
ble KAM-island in phase space which leads to chaotic 
scattering. A signature of this KAM-island is visible in 
representations of scattering wave functions in terms of 
Wigner functions. The avoided crossings of the Floquet 
quasi-energies of that system are studied in [171 ] . Sharp 
crossings of quasi-energies lead essentially only to a re- 
labeling of the states, whereas broad crossings, in which 
more than two states take part, completely alter the Flo- 
quet states and increase the high harmonic generation. 

Many other elementary systems have been studied as 
well: The so called tunneling diode, which consists of an 
oscillating quantum well between two static barriers, is 
investigated in and [ijj]. In this system the driv- 
ing creates additional sidebands due to the absorption 
or emission of oscillation quanta fko. The transmission 
through all sidebands exhibits a strong quenching at cer- 
tain parameters. The classical phase space of an infi- 
nite array of vertically oscillating potential wells is mixed 
chaotic at small energies and has regular spanning curves 
at high energies, see [l6j]. The motion of a particle in 
such a system resembles closely a random walk. Driven 
double- wells have been studied in [2(| and [2l|. The pe- 
riodic driving can be used in these systems to entirely 
suppress the tunneling between the two wells, if the driv- 
ing frequency is tuned to an exact crossing of the Floquet 
energies of the ground state doublet of the unperturbed 
double- well. 

Although several works deal with the tunneling 
through a laterally driven square potential, the classical 
dynamics in this system is largely unknown. The trans- 
mission of classical particles through a laterally driven 
Gaussian-shaped potential has already been calculated 
in However, this work does neither cover the whole 
range of possible parameters nor does it give any clue 
about the phase space structure, which is essential for 



the scattering process. The aim of this work is to close 
this gap and to provide a comprehensive survey of the 
periodically driven square barrier with oscillating posi- 
tion in the classical regime. We will analyze the entire 
phase space, where we find stable KAM islands, under- 
stand the underlying dynamics and the chaotic scattering 
process. We will also explore the full range of possible 
system parameters and make comparisons with quantum 
mechanical results for the transmission probability. 

This work is organized as follows: In section [TT] we in- 
troduce the system of the driven barrier and derive a 
mapping to describe the dynamics. In section IIIII we 
analyze the phase space where we find a stable island 
of quasi-periodic orbits. The position, size and parame- 
ter dependence of the stable island is studied in detail. 
The scattering dynamics, particularly the influence of the 
stable manifolds of the unstable periodic orbits, are the 
subject of section ITVl Due to the existence of the stable 
island, the scattering dynamics is chaotic and trapping in 
the sticky region of phase space is possible although the 
barrier is purely repulsive. Comparisons of the classical 
transmission probability with the quantum mechanical 
tunneling probability are made in subsection IIV Al Fi- 
nally, a summary is given in the last section fVl 



II. THE DRIVEN BARRIER AND ITS 
MAPPING 

Our classical system consists of a one-dimensional lat- 
erally oscillating potential of a finite and constant height 
Vq and width I, see Fig. [TJ The driving function is as- 
sumed to be harmonic, with a driving amplitude ao and 
frequency w 



V(x,t) = V e(- - |i-aocos(art)|) 
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FIG. 1: The ac-driven potential barrier. 
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Although the dynamics of the system is continuous, 
the forces acting on the particle are point-like and the 
particles move ballistically between collisions with either 
of the edges of the barrier. Therefore, it is sufficient to 
describe the dynamics in terms of a discrete mapping 
between collisions. To describe the particle-barrier in- 
teraction, we transform all coordinates into the frame of 
reference of the barrier, to the variables x and v. Even 
though this coordinate frame is accelerated and, there- 
fore, not an inertial frame, momentum and energy are 
conserved for the infinitesimally small time span of the 
interaction. 

x — > x — x — ao cos(wi) v — > v = v + oquj sin(cji) (2) 

When colliding with the barrier from the outside a par- 
ticle of mass to is transmitted into the barrier if its ki- 
netic energy, relative to the barrier, surpasses the barrier 
height Vq. 



E, 



171 ~2 
kin = ~r v n ^ *0 



(3) 



In this case the particle is decelerated due to energy con- 
servation 



Vn+l = \l — {Ekin — Vq). 
TO 



(4) 



If it is not transmitted, it is reflected and its new velocity 
becomes 



Vn+l = -V n - 



(5) 



Likewise, if a particle hits the barrier coming from its 
inside, it is always transmitted and accelerated 



v n +i 



— (Ekin +Vq). 



(0) 



Transforming these equations back to the laboratory 
frame yields the equations of motion 



v n +i = Vb(t n +i) + sign(v„ - v b (t n+1 )) 



(v n - V b (tn+l)) 2 ± —Vq, 

m 



(7) 



for transmission and 



v n+1 = 2v b - 



(8) 



if the particle is reflected, where v b (t) — — clquj sin(wt n +i) 
is the barrier's velocity at the time of the collision t n+ i 
and the sign ± depends on whether the particle is trans- 
mitted into the barrier (— ) or leaves the barrier (+). The 
time t n is mapped on the time t„+i of the next collision 
of the particle with one of the barrier's edges. Therefore 
i n +i is the smallest solution of 



Xb(t n +i) — x n + v n (t r 



(9) 



where x b can be either edge of the barrier, x b = o,q cos(wi) 
or x b = ciQCOs(ujt) + I. If equation ([9]) has no solutions 
for t n+ i > t n then the particle does not collide with 
the barrier again and escapes. This implicit equation 
can be solved only numerically. It is important to make 
sure that the numerically calculated collision time is the 
smallest solution of equation ([9]), because many effects, 
like stickiness in phase space (see sections IIIII and IIVI) , 
are susceptible to errors in the time mapping. 

The mapping shows that the five parameters of the sys- 
tem, the barrier's height Vq and thickness I, the driving 
frequency u>, the amplitude clq and the particle's mass 
to, can be reduced to just two effective parameters by an 
appropriate scaling transformation. Equations j7]) and 
(|8j then become: 

x -> x' = — t->t' = tuj v^v' = — (10) 

CLq UJCLq 



V 



n+l 



- sin(C+i) + si § n K + sin (C+i)) 



!<+sin« +1 )) 2 ± tt 



v 'n+i = -2sin(4+i) - v« 



(11) 

(12) 



The only parameter left in the mapping is where 
V u = Y a Q uj2 is the maximum kinetic energy a particle, 
which is at rest in the laboratory frame, can have in the 
barrier's frame of reference. The second parameter, 
is the barrier's thickness measured in units of the am- 
plitude and appears in the implicit equation ^ for the 
time mapping. In the following, we will scale the energy 
in units of V u and use the scaled coordinates in all cal- 
culations while keeping the same notation. The system 
is therefore completely described by the two parameters 
Vq and /. 



III. PHASE SPACE STRUCTURE 

We visualize the structures in phase space in terms of 
Poincare surface of sections by mapping all collisions of 
the particle with either side of the barrier to the Poincare 
section. This is equivalent to mapping all intersections 
of the trajectories in the 3-dimensional phase space with 
the two-dimensional manifold £1 defined by the barrier's 
motion 



CI = 




\t,ve$l 



where Xb(t) is the position of either of the barrier's edges. 
Since the driving function is assumed to be periodic, we 
can use the phase <p, defined as tot mod 2tt, instead of the 
time as a coordinate. The mapping naturally operates 
only on the manifold defined by CI, because it always 
maps a point in phase space on the point of the next 
collision with the barrier. The Poincare section is made 
unique by mapping only intersections with one chosen 
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edge of the barrier in a specified direction, i.e. v > Vb or 
v < vt, instead of all collisions. Thus the position and 
phase f are uniquely connected by the driving function 
and one of the two coordinates becomes redundant. In 
the following, we will discuss Poincare sections in which 
we plot the particle's velocity after a collision over the 
phase ip of the oscillating barrier. 

We covered the entire phase space with a fine grid of 
initial conditions to guarantee that all relevant structures 
are being shown in our Poincare sections. The resulting 
Poincare section is plotted in Fig. 
parameters are: Vq = 0.32 and I = 



2(a) 



0.4. 



The system's 
These parame- 
ters are typical for an experimental setup using semicon- 
ductor structures driven by external voltages or applied 
Laser fields, SSI. 




12 3 4 

phase of the collision 

(a) 




4.6 4.8 5 5.2 

phase of the collision 

(b) 



FIG. 2: (a): A typical Poincare section showing trapped par- 
ticles (b): Enlargement of island 4. 

The phase space of the ac-driven barrier has four sta- 
ble KAM-islands whose center is a stable periodic orbit 



of period four. This means that, through the driving, 
the repulsive potential can trap particles in a small part 
of phase space. This kind of dynamical trapping works 
only for a harmonic driving law or other similarly curved 
functions. The periodic orbit at the center of the island 
is stable because it lies symmetrically around the inflec- 
tion points of the harmonic driving function. A sawtooth 
shaped driving function, for example, does not create 
stable orbits. [39] An enlargement of the fourth island 



in Fig. 2(a) is plotted in Fig. 2(b) and shows the typi- 
cal structure of an elliptic fixed point surrounded by a 
stable island of quasi-periodic orbits and chains of sub- 
islands, which are the remnant of dissolved quasi-periodic 
orbits. The thin transition zone at the edge of the is- 
land is mixed chaotic and contains a fractal structure of 
sub-islands. The space outside the stable orbits contains 
only few points because this part of phase space is visited 
only by trajectories which leave the open system after a 
few collisions, whereas those on regular orbits stay in the 
scattering region (defined as the space that is covered 
by the barrier's oscillation) indefinitely. The Poincare 
section in Fig. |2(a) is not unique because it shows all 
collisions, with both edges of the barrier and in both di- 
rections. Structure 1 and 4 correspond to collisions with 
the left edge, structure 2 and 3 correspond to collisions 
with the right edge. The four regular islands are sym- 
metrical, the first and second structures are identical to 
the third and fourth with their phases increased by it 
and the sign of their velocity inverted. This reflects the 
symmetric properties of the driving function aocos(wi): 
cos(tt—lp) = cos(7r+<^) and cos(7r/2 — ip) = — cos(7r/2+iy9) 

The trajectory of the central periodic orbit of the sta- 
ble island, plotted in Fig. 3(a) is as following: Starting 
left of the barrier at (4) with a velocity of zero, the par- 
ticle is hit by the barrier at (1). Due to the relatively 
high negative velocity of the barrier at tpx, the particle 
is transmitted into the barrier and accelerated in nega- 
tive direction. This creates the stable fixed point of the 



island 1 in Fig. 2(a) and the corresponding stable island 
1 if we vary the initial conditions in the neighborhood of 
the fixed point. The barrier then overtakes the particle 
inside of it and the particle collides with the right edge of 
the barrier at (2). Since the velocity of the barrier at the 
points (1) and (2) is identical (vb(ipi) = Vb(f2)), the ve- 
locity of the particle becomes zero again. This collision 
corresponds to island 2 in Fig. 2(a) After the barrier 



has reached its minimum position and turning point at 
ip = 7T, it moves back in positive direction and hits the 
particle at (3). This collision accelerates the particle in 
positive direction and the particle moves with the bar- 
rier until it is overtaken by it at (4), where the velocity 
of the particle becomes zero again and the periodic cycle 
starts again. Due to the symmetry of the driving func- 
tion (cos((p)), the collisions (1) and (2) are symmetrical 
to (3) and (4). 

The trajectories of the quasi-periodic orbits can be un- 
derstood as a perturbation of the periodic orbit described 
above. A typical trajectory of such an orbit is plotted in 
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FIG. 3: Motion of a trapped particles. The black line rep- 
resents the particle, the two gray lines represent the edges of 
the barrier, (a) is the trajectory associated with the central 
periodic orbit of the stable island, (b) is a typical trajectory 
of a quasi-periodic orbit. The numbers in (a) correspond to 
the island numbers in Fig. 2(a) 



Fig. 3(b) and shows both modes of the motion: the pe- 



riodic hopping between the four islands inherited form 
the stable periodic orbit and an overlaid harmonic oscil- 
lation. This sinusoidal motion is represented as closed 
orbit in the Poincare section. The frequency of this har- 
monic oscillation, measured by a Fourier transformation, 
is equal to the frequency by which the particles rotate on 
the quasi-periodic orbit around the elliptic fixed points 
in the Poincare section. 

It is possible to calculate the position of the central 
elliptic fixed points, here we will do it for island num- 
ber 1. As Fig. 3(a) shows, the velocity of the periodic 



orbit is zero while the particle is outside of the barrier 
in structure 2 and 4 (t>o = «2 = 0). Due to the sym- 
metry of the system, the collision points are symmetric 
around \ in their phases ip± and tp2 (<f2 — 7t — ipi) and 
around the equilibrium position in their positions x\ and 
X2- Therefore, the velocity in structure 1 has to be: 



Ax 
At 



(f -</>!) 



x b (ipi) = cos(^i) (13) 



where Xb is the position of the left edge of the barrier. 
This velocity can be calculated from equation 



vi = v b {ipi) + sign (v Q - v b {Lpi)) y (v - v b {ipi)) 2 - V 

The initial velocity vo is zero and v b = — sm(<p). The 
result is an implicit equation for tpy. 



2 cos(y>i) — I \ " f 2 sm(ipi)(n — 2y?i) 



7T — 2lfi 



2 cos(yi) — I 



- 1 



Vb 



(15) 



This implicit equation determines the phase of the 
first elliptic point. For values of I e [0,2], the equa- 
tion f(ip) = Vo has two solutions in the interval [0, ?]. 
The physically relevant solution lies to the right of the 
function's maximum. For Vq > f m ax, where f m ax is the 
maximum value of / at a chosen value of I, the implicit 
equation (|15p has no solutions. This means that the bar- 
rier's potential is too strong to allow for trapped orbits. 
The function has no positive values for I > 2 and the sec- 
ond root is greater than ^ for I = 0. The elliptic orbits 
disappear in both cases. Figure [4] shows the maximal 
value of f(<p) as a function of I. Only pairs of parameters 
(Vb, I) below the curve allow for trapping of particles. For 
sets of parameters (Vb,Z) above the curve in Fig. Q] the 
phase space contains no bound orbits. To be precise, the 
central periodic orbit does not just become unstable for 
other parameters. The periodic orbit, the surrounding 
elliptic island and all unstable periodic orbits cease to 
exist! 

The maximal values of Vb = 1 and I = 2 can be easily 
understood physically: is the maximum kinetic energy 
a particle can have in the barrier's frame of reference if 
it is at rest in the laboratory frame. If Vb is greater 
than Vu,, then particles with velocity zero will never be 
transmitted into the barrier. Since all stable orbits cross 
the v = axis, this would destroy all of them. For I > 2, 
t>i in equation (|13p would have to be positive, which is 
physically impossible (see discussion of the periodic orbit, 
Fig. [3(a) |. 



The phases of the other fixed points can be calculated 
from the phase and position position of the fixed point of 
island 1 by using the symmetry properties of the driving 
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FIG. 4: Maximum barrier potential allowing for periodic or- 
bits as a function of the barrier's thickness. Only pairs of 
parameters below this curve lead to stable orbits. 



function: 

lfi 3 = 7T + Ifi 

9?4 = 2tt — Ifi 

The shape and size of the stable islands in the Poincare 
sections change with varying parameters. Fig. 2(b) 



is 

typical for the phase space of this system: The central 
fixed point and the island of quasi-periodic orbits are sur- 
rounded by a chaotic layer with a fractal structure of sta- 
ble and unstable periodic orbits. Additionally, there exist 
one or more sets of large and distinguished sub-islands. 
These can be inside the main island, as the five sub- 
islands in Fig. |2(b)| or outside of it. These sub-islands 
discern themselves from other KAM substructures in sev- 
eral ways. They are not only much larger than other sub- 
structures, but their creation and destruction follows a 
simple pattern as the parameters are changed. Keeping 
I constant while decreasing Vo, the sub- islands of period 
n move towards the outer edge of the stable island. As 
Vo is decreased further, the sub-islands cross the out- 
ermost quasi-periodic orbit of the main island, forming 
separate sub-islands, and ultimately dissolve. Simulta- 
neously, a new set of sub-islands of period n + 1 forms 
at the center of the stable region and begins to move to 
the outside. Islands with an even period appear as pairs, 
therefore the sequence of the number of large sub- islands 
is 4,3,8,5,12,7,..., see Fig. [5] The sequence starts at 
the maximal value of Vb that allows for bound orbits as 
plotted in Fig. [5] In the limit of very small Vb the period 
n diverges and the sub-islands form almost a continuum 
that can not be resolved numerically. Between period 
n = 4 and n = 3 the KAM island takes a triangular form 
and its size goes to zero. The creation and destruction 



of the sub-islands follows the same sequence for all I as 
Vo is varied, albeit on a different scale. A variation of 
the barrier width I at a constant Vq leads to the same 
sequence as well. This suggests that the qualitative be- 
havior of the phase space structure can be described by 
just one effective nonlinearity parameter, consisting of a 
combination of I and Vo. 

This sequence of sub-islands is typical for nonlinear 
systems and has been studied in detail in the standard 
map, see e.g. [U, [25|. The creation and destruction of 
stable periodic orbits, or rather quasi-periodic orbits in 
general, is tied closely to number theory. The KAM tori 
can be characterized by their winding number. It was 
conjectured [26| that the last KAM tori to be destroyed 
when the nonlinearity is increased are those with rotation 
numbers equal to noble numbers, which can be written 
as continued fractions [4(j. Before such a torus with ro- 
tation number R is destroyed, all the periodic orbits with 
rotation number equal to a truncation of the noble num- 
ber R become unstable. Thus, the winding numbers of 
the large sub-islands can be interpreted as the trunca- 
tions of the simplest first-order noble numbers. 

These sub-islands also play a significant role in the 
scattering process. As we will demonstrate, the flow into 
and out of the border zone of the stable island is dom- 
inated by the stable and unstable manifolds of the out- 
ermost large sub-islands, even for parameters for which 
the island has entirely dissolved. 

The area covered by the elliptic island changes with 
the parameters as well. We analyzed numerically the 
surface covered by the stable island as a function of the 
parameters by dividing the phase space of the Poincare 
sections into a fine grid of 10 6 small squares. All squares 
containing data points of the stable islands count as part 
of the surface. We checked this method for its stability 
by doubling the number of data points and comparing 
the resulting surfaces. The corresponding area is shown 
in Fig. [5] as a function of the barrier's height Vo for dif- 
ferent values of the barrier's width I. The peaks of the 
area in Fig. [6] are well understood: The size of the stable 
island is defined by two competing effects. As the poten- 
tial hight Vo is decreased, the quasi-periodic orbits move 
outwards, which increases the island's size. At the same 
time, outer quasi-periodic orbits are destroyed. Since tori 
with noble rotation number are destroyed very late, the 
covered surface reaches a maximum whenever the outer- 
most torus has a noble rotation number. At this point, 
other tori inside this noble torus have already been de- 
stroyed and the outermost curve separates a chaotic sea 
inside of the island from the outer chaotic area. When 
this outermost noble torus is destroyed, the two chaotic 
parts become connected and the size of the island reaches 
a sharp minimum. This kind of behavior is typical for 
nonlinear systems, see [2ol [27|. 

To gain further insight into the properties of the tran- 
sition zone around the stable island, we use a method 
developed in [28| and [2£| to locate the unstable periodic 
orbits (UPO) of this system. Given a N-dimensional dis- 
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FIG. 5: Overview of the parameter dependence of the KAM- 
island. The parameters are I = 0.1 and (a): Vo = 0.895 (b): 
Vo = 0.867 (c): Vo = 0.86 (d): V = 0.76 (e): V = 0.7 (f): 
Vo = 0.553 (g): Vo = 0.5 (h): V = 0.05 



Crete chaotic dynamical system U defined by 
U : f l+1 = f(n) 



(16) 



a linear transformation is used to construct a new system 
Sk defined as 



Sk : r i+1 =fi + \C k (f p (ri) - r*) 



(17) 



where the matrix Cfe is orthogonal, A is a small factor 
and p is the period of the fixed point we want to stabi- 
lize. Evidently, Sk and f p have the same fixed points. It 
can be shown that for every unstable fixed point there ex- 
ists a suitable orthogonal transformation matrix Cfc and 
a small factor A such that this unstable fixed point is 




--1=0.1 
1=0.4 
- 1=0.8 
— 1=1 .2 



0.8 



FIG. 6: Phase space volume of the stable island as a function 
of Vo for different I. 



stable under Sk- The factor A has to be small enough so 
that the eigenvalues of the matrix 1 + ACfc(TV — 1) have 
absolute values smaller than 1, where TV is the stability 
matrix of the system U. The matrices Cfc correspond 
to reflections and rotations along the coordinate axes, 
thus all entries of Cfc are Cy £ {0, ±1} and each row 
and column contains only one element which is different 
form zero. There exists a total number of N\2 N of such 
matrices that will, in general, stabilize different types 
of unstable periodic orbits. However, it can be shown 
that a much smaller number of Cfc is sufficient to find 
all periodic orbits because some of the linear transfor- 
mations are redundant. Only three matrices are needed 
in a two-dimensional system, see 30] . The advantage of 
this method over more conventional approaches, such as 
Newton-Raphson, is the global convergence. Even ini- 
tial conditions far away from a fixed point eventually 
converge if the matrix Cfc and the factor A are chosen 
appropriately. 

In the system of the driven barrier we search for pe- 
riodic orbits in the unique Poincare sections such as 
This way each point in the (v, (^)-plane is 



Fig. 2(b) 



uniquely connected to a point in the (x, v, t) phase space. 
To find the periodic orbits we cover the (v, (^)-plane with 
a grid of 10 initial conditions and iterate the transformed 
mapping with all three matrices Cfc and for periods of one 
to 21 with respect to the Poincare section. The parame- 
ter A is chosen between 5 • 10~ 3 and 10~ 5 for higher peri- 
ods. A typical result is plotted in Fig. [7] which shows the 
unstable periodic points as black crosses, the stable peri- 
odic orbits as black dots and the corresponding Poincare 
section in gray. The stability was derived according to 
the eigenvalues of the monodromy matrix. The chaotic 
layer around the elliptic island contains many families 
of periodic orbits which form the skeleton of the frac- 
tal structure. The number of unstable fixed points rises 
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exponentially with their period. n „ 
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FIG. 7: Plot of the Poincare section and the periodic orbits 
of period one to 21. Stable orbits are shown as dots, unstable 
orbits as crosses. 

The unstable periodic orbits, or rather their stable 
manifolds, play an important role for the scattering pro- 
cess. This is because, although the stable island is itself 
not directly accessible by initial conditions starting out- 
side of the interaction area, the stable manifolds, also 
called stable asymptotic curves, reach far into the part 
of phase space which is accessible from the outside. In 
order to calculate the flow of an unstable periodic or- 
bit we take a small initial segment of length 10~ 8 along 
the eigenvectors of the monodromy matrix at the UPO 
and iterate an ensemble of 10 6 initial conditions in this 
segment forward in time along the unstable asymptotic 
curves or backward in time along the stable asymptotic 
curves. It is in this system not possible to follow the man- 
ifolds of a specific UPO, because the manifolds of differ- 
ent UPOs cross each other at heteroclinic intersections. 
Thus the flow becomes chaotic and our simulations pro- 
duce a global picture of the flow of the system of UPOs. 
The stable and unstable manifolds, plotted in Fig. 8(a) 
and Fig. |8(b)| reach far out of the border region into the 
chaotic sea in what looks in these plots as spiral arms, 
which consist of an infinite number of stable or unstable 
curves. These outer manifolds belong to periodic orbits 
with a period of 3 or multiples of 3. 

We calculated the asymptotic curves of the system for 
the whole range of parameters and found that they are 
closely tied to the large primary sub-islands around the 
island. The flow into and out of the stable island is al- 
ways dominated by the outermost family of sub-islands, 
even if this family is completely unstable. Therefore the 
number of spiral arms which are formed by the mani- 
folds follows the same sequence as the sub-islands when 
the parameters are changed. The unstable periodic or- 
bits and with them their asymptotic curves are present 
only for parameters which allow for the existence of the 




(b) 



FIG. 8: Stable (a) and unstable (b) manifolds of the UPOs 
of the system. 



stable island in phase space, see Fig. 2J 



IV. SCATTERING DYNAMICS 

The phase space structure presented in section UTT1 has 
profound effects on scattering processes. Due to the ex- 
istence of a KAM-island in phase space, the ac-driven 
barrier is a chaotic scatterer. (See [31| for a definition.) 
To simulate the scattering process we place an ensemble 
of particles with an uniformly distributed velocity Vi n far 
outside of the scattering region. The initial phase tpo (or 
time) is distributed in such a way that the phase of the 
first collision tpi with the barrier is uniformly distributed 
in [0, 27r]. It makes sense to use the initial velocity Vi„ 
and the phase of the first collision tpi as parameters of 
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the scattering because this allows an easy comparison to 
the Poincare sections of section IIII1 where we use simi- 
lar coordinates. The disadvantage of these coordinates 
is that not all combinations of Vi n and tpi are accessi- 
ble from the outside of the scattering region. Thus, the 
following plots have an area marked as inaccessible. 

As scattering functions we examine in detail the dwell- 
time, defined as the time the particles spend in the scat- 
tering region, and the change of the velocity of the parti- 
cles, \vi n \ — \v out \. We also determine the particles' num- 
ber of collisions and whether the particles are transmit- 
ted or reflected. The velocity change is plotted in Fig. [5J 
for Vq = 0.32 and I = 0.4 in a gray-scale plot. Dark 
surface colors represent acceleration, bright areas stand 
for deceleration. In the following, we analyze the scat- 
tering process for this representative set of parameters. 
It is possible to discern different well separated regions 
in Fig. [5] in which the scattering function is smooth. In 
other regions, most prominently in a wedge-shaped part 
in the center of Fig. the scattering function has unre- 
solved parts. The scattering process is chaotic in these 
regions. 




2 4 
phase of the first collision 



FIG. 9: Velocity change |«j n | — \v out \ of the particle as a 
function of velocity Vi n and phase of the first collision <pi in 
gray-scale. 

We discuss the case of regular scattering first: The ve- 
locity of fast particles is hardly changed at all, because 
fast particles simply traverse the barrier in a very short 
time. According to eq. j7]), the transmission through 
the barrier is elastic if the velocity of the barrier at the 
collision with the left edge of the barrier is equal to 
the barrier velocity at the collision with the right edge, 
v b(<Pi) = v b( i P2), which is approximately the case for fast 
particles. The scattering of slow particles is, in general, 
inelastic, particularly if a particle is reflected by the bar- 

and |10(bj| show trajectories 



Figs. 10(a 



rier, see eq. 

typical for acceleration and deceleration, respectively. 

The regular regions visible in Fig. [5J correspond to a 
constant number of collisions, as plotted in Fig. [TTJ The 
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FIG. 10: Typical trajectories for acceleration (a) and decel- 
eration (b) at (vi n = 0.4,ipi = 2.3) and (w,: n = 1.5,y>i = 4.5) 



edges of the regular regions correspond to a change in 
the number of collisions, not necessarily by one, which 
naturally leads to a sudden change in the other scattering 
functions. The number of collisions is small, between one 
to four, for all regular scattering trajectories. 
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FIG. 11: Number of collisions of the particle as a function 
of velocity Vi n and phase (pi represented as gray-scale. The 
collision number is also printed in the plot. 

In the chaotic parts of the scattering function in Fig. [51 
the dynamics is infinitely sensitive on the initial condi- 
tions. Therefore it is not possible to resolve the scattering 
function in these areas completely. These unresolvable 
points are singularities. We have tested this by succes- 
sively magnifying the irregular parts of Fig. [9] up to the 
numerical limits. Fig. 1121 shows such an enlargement by 
a factor of 10 13 . The singularities of the scattering func- 
tion, visible as unresolved parts in Fig. [T2] and [5J form 
a fractal set. Actually, the chaotic parts of the scatter- 
ing function are predominant on smaller scales, i.e. the 
smooth parts, still visible in Fig. [T21 become rare. This 
is typical for non- hyperbolic chaotic scattering ([IH). 

The origin of the chaotic scattering is the stickiness of 
particles to the KAM-island in phase space. This can be 
seen from a logarithmic plot of the dwelltime in Fig. 1131 
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FIG. 12: Enlargement of part of Fig. H by a factor of 10 13 . 



The dwelltime diverges in exactly the regions in which 
the scattering function has singularities. The number of 
collisions diverges as well in these chaotic regions. (This 
can not really be seen in Fig. 111! because the color map 
is capped at n = 10. But the collision number reaches 
n = 10 4 and more in the white parts of the plot.) Trajec- 
tories starting on initial conditions in the chaotic parts of 
the scattering function become sticky and trace the out- 
ermost quasi-periodic orbits of the stable island for arbi- 
trary long times. The number of collisions per time unit 
is therefore constant for all sticky trajectories, namely 
four collisions per period T of the driving function. 
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FIG. 13: Logarithm of the dwelltime of the particle as a func- 
tion of velocity Vi„ and phase ipi represented as gray-scale. 

It is known that in most systems which exhibit chaotic 
scattering, singularities in the scattering function have 



a divergent dwelltime. (|3i() This can be easily under- 
stood since the scattering function is infinitely sensitive 
to perturbations of the initial conditions leading to sin- 
gularities which can only be the case if the interaction 
time in continuous systems or the number of interactions 
in discrete systems between target and particle diverges 
as well. 

The dwelltimes of sticky particles have a typical prob- 
ability distribution. Fig. [14] shows the distribution of 
dwelltimes in the sticky regions of Fig. [T2] for a total of 
more than 10 10 random initial conditions. The dwell- 
time distribution can be approximated by a power law 
P(td) ~ td- We find an exponent of 7 = —2.5 from a fit 
to Fig. [H between t d = 10 1 and t d = 10 8 . 
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FIG. 14: The distribution of the sticking times, using a log- 
arithmic adjusted bin size. By fitting a power law we found 
an exponent of 7 = —2.5. 

The origin of the stickiness in this system are the sta- 
ble manifolds of the unstable periodic orbits described in 
section IIIII Although the KAM- island itself is inaccessi- 
ble by scattering trajectories, the stable manifolds reach 
out of the inaccessible area and intersect the set of scat- 
tering initial conditions. These intersections are identical 
to the singularities of the scattering function and have an 
infinite time delay. The fractal structure of the scatter- 
ing function in the chaotic regions is just the structure of 
the stable manifolds, which arc formed by the unstable 
periodic orbits in the KAM-island. 

Another way of characterizing a chaotic scattering pro- 
cess is to calculate the uncertainty dimension of the 
chaotic part of the scattering function, which is a proxy 
for the fractal dimension. The uncertainty dimension de- 
scribes the scaling with the resolution of the proportion 
of singular points to regular points in a two-dimensional 
scattering function. As predicted in [32| for nonhyper- 
bolic chaotic scattering systems, such as the ac-driven 
barrier, the uncertainty dimension is exactly one for all 
parameters that allow for the existence of the stable 
KAM island. 
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The scattering function has two other singularities, 
which are not caused by stickiness to regular structures. 
Such singularities are isolated and do not lead to chaotic 
scattering. Nevertheless, the scattering function is in- 
finitely sensitive to changes in the initial conditions at 
these singularities. One of them is located at Vi n = 1 and 
(pi = 1.5tt in Fig. [5] and is called "whispering gallery" in 
static systems. The dwelltime, see Fig. [T3J is quite small 
at this point, td = T/2, whereas the number of collisions 
diverges, see Fig. QT] A particle on this trajectory hits 
the barrier at tpi = 1.5tt with a velocity slightly larger 
than 1, Vi n = 1 + e. The barrier has its maximum ve- 
locity of Uft = 1 at this phase. The particle is reflected 
and the new velocity becomes Vi — 2v b — v- m = 1 — e 
according to eq. ((8|). Because the barrier decelerates 
after cp = l.bn the particle collides with the barrier 
again after a very short time at if2 ~ 1.57T + V6e. 
The particle is decelerated by this collision further to 
v 2 » 2v b -vi = 2(1 - 3e) - (1 - e) = 1 - 5e. This se- 
quence of successive collisions continues until the barrier 
reaches the turning point at <p = 2ir from where, due 
to the symmetry of the driving function, the process is 
inverted and the particle is accelerated by successive colli- 
sions. Since the particle has the last collision at <p = O.bn 
the dwelltime is only td = T/2. However, for sufficiently 
small e the particle can have an arbitrarily large number 
of collisions, because the time passed between the colli- 
sion is proportional to e. A typical trajectory is plotted 



in Fig. 15(a) 



The other singularity is a so called low velocity peak, 
[331 ]. It only appears if the width of the barrier is larger 
than twice the oscillation amplitude, I > 2, and if the 
potential height is greater than the effective potential, 
Vq > 1. For these parameters, a particle can hit the 
barrier at the extremal position in tp = at a velocity 
of exactly vq = \/Vq + e. Such a particle is transmitted 
into the barrier and is decelerated to the small velocity 
Vi = y/e. Due to the large width of the barrier, the 
dwelltime of the particle becomes approximately td m 
and thus diverges. A typical trajectory of a low velocity 
peak is plotted in Fig. |15(b)| 
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FIG. 15: Singularities of the scattering function not associ- 
ated with chaotic scattering, (a) is the trajectory of a whis- 
pering gallery orbit, (b) is a low velocity peak. 

As the parameters of the system are changed, the regu- 
lar parts of the scattering function (Fig. [5]) are deformed 



and shifted to different initial velocities, but remain qual- 
itatively similar. The shape and position of the chaotic 
parts stay qualitatively similar as well, whereas their ex- 
act structure depends critically on the parameters. This 
is because the chaotic regions are created by the stable 
manifolds of the UPOs. The shape of the stable mani- 
folds is directly connected to the primary sub-islands and 
therefore follows the sequence of creation and destruction 
described in section |TTT] as the parameters are changed. 
Since the stable island and the surrounding unstable pe- 
riodic orbits are the cause of the chaotic scattering, the 
scattering on the ac-driven barrier becomes regular for 
parameters that don't allow for elliptic orbits in phase 
space (see Fig. U]). The singularity corresponding to the 
whispering gallery is independent of the parameters of 
the system, it exists in all harmonically laterally driven 
systems and it is not a unique property of the driven 
barrier. The low velocity peaks require a large barrier 
height Vb > 1 and width I > 2. Consequently, there are 
no parameters for which the low velocity peaks and the 
chaotic scattering due to the KAM-island coexist. It was 
shown in Refs. [H, [H[ that low velocity peaks can lead 
to a new form of scattering dynamics, called dilute chaos, 
which requires the existence of UPOs and therefore does 
not appear in this system. 



Application and comparison to the quantum 
behavior 



We can use our knowledge of the system to make com- 
parisons with established results for the dynamics of the 
driven barrier in the quantum regime. The tunneling 
through a periodically driven square potential barrier has 
been analyzed in [ll[. (See also Q.) It was found that 
the transmission coefficient as a function of the parti- 
cle energy has resonances below the minimal tunneling 
energy of the static system for high driving frequencies. 
This can be explained by resonant tunneling into semi- 
stable bound states of an effective time-averaged poten- 
tial, which has a double-barrier structure. We now want 
to compare the results in the quantum regime with our 
classical simulations. The parameters chosen in [ll| are 
a = 200, I = 80, m = 0.1 and V = 0.0147. The sys- 
tem is studied for three values of the driving frequency, 
u) = 0, u) = 3 • 10~ 4 and u> = 3 • 10~ 2 . The corresponding 
effective parameters are I — 0.4 and Vb — > 00, Vb = 81.6 
and Vb = 0.0082. To imitate the transmission of a quan- 
tum wave packet with a classical simulation we use an 
ensemble of particles that is distributed as a minimal un- 
certainty Gaussian in position and momentum space, i.e. 
PxCp = h- The results of our classical simulation are 
shown in Fig. [16] 

Surprisingly, the classical simulation coincides to an 
amazing degree with the results of the quantum mechan- 
ical analysis in the high frequency limit (Fig. 16(c) I. The 



transmission in the classical system shows the same res- 
onant behavior as in the quantum regime. (It should be 
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FIG. 16: Transmission function for different driving frequen- 
cies of the classical system, (a): lo — 0, (b): u = 0.0003, (c): 
lo = 0.03 



mentioned that a fine tuning of the initial conditions is 
necessary to reproduce the resonances of [ll| exactly.) 
Since the model of the effective time-averaged potential 
V e f f fails in the classical regime, this result is surprising. 
The scattering of classical particles off the static potential 
V e ff would simply reproduce a step function, since the 
classical mechanics just does not allow any tunneling into 
resonant states. It is important to note that the effective 
potential is generally ill suited to describe the dynamics 
of trapped particles in the classical regime. Although the 
effective potential could in principle be used to explain 
the existence of trapped particles for high frequencies, we 
also find trapped particles for low frequencies, i.e. when 
the driving frequency and the oscillation frequency of the 
trapped particles are of similar order of magnitude. The 
trapping is not caused by an effective potential for high 
driving frequencies but by a synchronization of the mo- 
tion of the particle and the barrier. 

The origin of the resonances of the transmission func- 
tion in the classical regime is very different from the 
quantum regime. The transmission function of a sin- 
gle particle in Fig. 17(a) does not show any resonances 
as a function of the initial energy, at least not of the 



form seen in Fig. 16(c) (It makes sense to compare the 



transmission function of a single particle Fig. 
Fig. ~ 



17(a) with 
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FIG. 17: Transmission function of a single particle as a func- 
tion of the initial velocity and (a): the phase of the first col- 
lision ifii (b): the initial phase tpo- 



ble changes the qualitative behavior very little.) The 
reason is that the phase of the first collision ipi, used 
in Fig. 17(a) as coordinate, is not an appropriate coordi- 



nate to describe this scattering process. In the simulation 
leading to the transmission resonances in Fig. 16(c) the 



16(c) because using a narrowly distributed ensem- 



initial phase ipo is kept constant whereas (pi oscillates 
wildly as a function of the energy E. When the initial 
velocity Vi n of a particle starting at a distance of Xo from 
the scattering region is changed by a small amount Av, 
the time when the particle enters the scattering region 
is changed by At = -%g-Av. In all these simulations, the 
initial velocity v in is very slow, therefore the variation 
of the initial energy E leads to a fast oscillation of the 
collision phase tp±. When we plot the transmission as a 
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function of the initial phase in Fig. |17(b)| the resonances 
of Fig. 16(c) become visible. 



Whether an incoming particle is transmitted or re- 
flected depends, for small particle energies, mostly on 
whether the barrier is approaching the particle or re- 
ceding from it at the time when the collision occurs. 
A variation of the energy for a constant collision phase 
(pi has little impact on the transmission probability, see 
Fig. 17(a) Therefore, these resonances are not really res- 



onances of the energy of the particle, as in the quantum 
regime. They are produced simply by the propagation to 
the barrier, where the different collision phases lead to 
peaks in the transmission probability. 

It would be very interesting to analyze the similarity 
between the classical and the quantum case in the frame- 
work of semiclassical physics. Doing so would require ap- 
plying the Gutzwiller formula to the periodic orbits of the 
system to derive a semiclassical propagator for the sys- 
tem, see (36[. However, such a study is beyond the scope 
of this work. It should be noted that in related systems 
such as an oscillating hard wall pot ential, semiclassical 
approximations did not work [37l. |38|| . 



V. SUMMARY AND OUTLOOK 

The aim of this work is to analyze and understand in 
detail the classical dynamics of the ac-driven barrier for 
the full range of parameters. Although the potential is 
repulsive, the system exhibits a dynamical trapping pro- 
cess which is associated with an island of stability. This 
trapping process can be understood as a synchronization 
process between particle and barrier, which depends on 
the curvature of the driving law. The stable KAM-island 
of quasi-periodic orbits in phase space leads to topologi- 
cal chaos. The central periodic orbit, and with it all sta- 
ble and unstable periodic orbits, exist only for a limited 



range of parameters. We determined these parameter 
ranges, calculated the position of the period four orbit, 
the size of the elliptic island and its shape as a function 
of the parameters. The transition zone around the stable 
island contains an infinite set of unstable periodic orbits, 
the stable manifolds of which reach far away from the 
stable island. These stable manifolds make the system a 
chaotic scatterer. Initial conditions starting on the stable 
manifolds are singularities and have a divergent dwell- 
time and collision number. The singularities form a frac- 
tal set with an uncertainty dimension of one. The system 
possesses two additional types of singularities, the whis- 
pering gallery and a low velocity peak. These are isolated 
singularities and are not connected to the KAM-structure 
in phase space. The transmission function of a suitably 
prepared ensemble yields results which are very similar to 
tunneling resonances in the quantum mechanical regime. 
However, the origin of these resonances is very different 
in the classical regime and this sheds a new light on the 
high frequency behavior of the driven barrier. 

The results of this work all depend on the existence of a 
dynamical trapping process. The stable orbits which we 
discovered rely on the curvature of the harmonic driving 
law. When we use a sawtooth shaped driving law, we 
find no such stable orbits. It can be assumed that other 
suitably curved driving functions lead to bounded motion 
as well. Likewise, our results do not depend on the exact 
shape of the barrier itself. 
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